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A FUNCTIONAL CLT FOR THE OCCUPATION TIME OF A 
STATE-DEPENDENT BRANCHING RANDOM WALK 

By Matthias Birkner and Iljana Zahle^ 
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and Universitdt Erlangen-Niirnberg 

We show that the centred occupation time process of the origin 
of a system of critical binary branching random walks in dimension 
d > 3, started off either from a Poisson field or in equilibrium, when 
suitably normalized, converges to a Brownian motion in d > 4. In 
d — 3, the limit process is a fractional Brownian motion with Hurst 
parameter 3/4 when starting in equilibrium, and a related Gaussian 
process when starting from a Poisson field. For (dependent) branch- 
ing random walks with state dependent branching rate we obtain 
convergence in f.d.d. to the same limit process, and for d — 3 also a 
functional limit theorem. 



1. Introduction and main result. We study the fluctuation behavior of 
the occupation time in a single point of a system of critical binary branching 
random walks with a state-dependent branching rate (BRW). BRW consists 
of particles which move independently on Z*^ in continuous time according 
to a given random walk kernel a. The branching rate at a site depends on 
the number of particles there: if there are k individuals at x, a branching 
event in x occurs at rate a{k). The particle which is chosen to branch then 
leaves either two or zero offspring at its current location, each possibility 
occurring with probability 1/2. The classical case of independent branching 
with constant branching rate p corresponds to cr{k) = pk. We assume further 
on that the branching rate function a is not = and Lipschitz: 

(1.1) \a{m) — a{n)\ < c\m — n\, 
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which especiahy imphes 

(1.2) a{k) < C2k for some < C2 < oo. 

We denote by Ct(x) the number of particles present at location x at time 
t. We assume that the transition rate matrix a{x,y) = a{0,y — x) governing 
the individual motion of particles is symmetric, irreducible and has finite 
second moments, which implies 

(1.3) {Qij)i,j=i,...,d = [ ^ a(0,x)xjXj is finite and invertible, 

VxeZ'^ / i,j=l,..;d 

where x = (xj)j=i^...^c;. We have J^x'^i^^^)^ = by symmetry, and we can 
assume without loss of generality that a is stochastic, that is, ^^(0) x) = 1. 
Denote the continuous time transition probabilities by at{x,y). 

It is well known that the independent BRW in d<2, starting from any 
initial condition with bounded local density, suffers local extinction, that is, 
£,t{x) — > in probability as t — > oo for any x S Z"^. For translation invariant, 
ergodic initial distributions with finite intensity, this can be found in, for 
example, [17]. It can be extended by the comparison argument in the proof 
of (5.1) in [13]. 

On the other hand, in d>3, there exists a one-parameter family of ex- 
tremal invariant probability measures A^, "!? > 0, parametrized by the ex- 
pected density: / ^{x)A^{d£,) = t). Each is shift-invariant, and {^t{x) : x £ 
U^^t > 0} under A^ is ergodic with respect to space- and time-shifts. See 
references below Theorem 2.3. 

Let us denote the distribution of a Poisson field on Z*^ with homogeneous 
intensity 'd by 7{{'d), that is, under 'H{'&), the random variables ^(x), x £ Z'^, 
are i.i.d. Poisson('i?). If -C(^o) = T~(-{'&), we have C{(,t) — > A^ weakly as t — > oo. 

First we present the main result for the case of independent branching 
[i.e., a{k) = pk]. Let >C(Co) G {W('!?), A^}. By ergodicity, the occupation time 
of any point x G Z'^ satisfies 

1 

— / Ct{x) — > -i? almost surely as T ^ oo. 
T Jo 

Thus, a natural question concerns the random fiuctuations of the occupation 
time around its asymptotic limit. This is the content of our main result: 



Theorem 1.1 (Independent branching). 1. If {£,s)s>o started in the 
(unique extremal) equilibrium distribution with intensity > 0, then the 
processes 
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converge in distribution as elements o/C([0, oo), M) (denoting the set of con- 
tinuous functions on [0, oo) with values in M) toward a Brownian motion 
in d>4: and to a fractional Brownian motion with Hurst parameter 3/4: in 
d = 3 as N ^ oo, where the norming is given by 

(t^/\ d = 3, 

(1.5) hd{t) = \ ^/m^t, d = A, 

The covariance of the limiting process X is given by 
Coy{X,,Xt) 
V2 



3^3/,(detQ)-V2^p[t3/2 + ,3/2_|^_^|3/2j^ 
(27r)-2(detg)-i/2i9px (sAt), 



duau{0,0) + p duuau{0,0) 



X (s At), 



d = 3, 
d = 4, 
d>5. 



2. The same conclusions hold if C{^q) = 7i{'d), and d > 4. In the case 
C{^q) = 'H{'d) and d = 3, the processes X^ converge toward a Gaussian pro- 
cess X with covariance given by 

Cow{Xs,Xt) = |^(detQ)^i/2^p 

(1.6) 



t3/2 + s3/2_i|t 



1 3/2 



-it + s 



,3/2 



The normalizations h^ are dictated by the requirement of a nontrivial 
covariance function for the hmit process, and this in turn is determined 
by the decay properties of the transition probabihties of the underlying 
random walk a; see the calculations in Section 3.3. Note that with /) = 0, 
BRW becomes a system of independent random walks, and has the family 
7i{'d), "i? > 0, of shift-invariant extremal equilibria. In the situation p = 0, we 
see from Theorem 1.1 that under the normalizations used in Theorem 1.1 
the limit process X is trivial in d < 4 and a Brownian motion in d > 5. This 
is in keeping with a "metatheorem" that the introduction of branching shifts 
"critical dimensions" by 2: In a system of independent random walks, the 
occupation time requires normalization by t^^^ in d= 1 , \/ 1 log t in d = 2 and 
^/i in (i > 3 in order to obtain a nontrivial limit (see [12]). 

While for nonbranching random walks, the nonclassical norming is due 
to recurrence properties of the individual particles, the behavior in our case 
is governed by the recurrence properties of families: The equilibrium of a 
BRW can be decomposed into a Poisson system of "clans" of particles with 
a common ancestor (see, e.g., [28]), and such a clan will visit the origin in- 
finitely often if and only if d < 4 (see [26], Theorem 1). This allows in the case 
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of independent branching, at least on a heuristic level, also to understand 
the different normings. Substituting the probabilistic representation of the 
Palm distribution at at time T from [26], Proposition 1 for "a typical clan 
which contributes to the occupation time up to T," we see that the expected 
contribution per clan is 

(1.7) r dt [' dsaT+t-2s{0,0), 

Jo J A 

where ^ = when C{(,q) = Ti-i'd) and A = —oo when >C(^o) = ^i?- In d = 3 
this grows like VT, hence, we expect of the order T/^/T = y/T clans to con- 
tribute. So due to independence of clans, the fluctuations should indeed be 
of order \/T(v^)"^'^^ = T^/^. In d = 4, this grows like logT, so T/logT clans 
should contribute, suggesting fluctuations of the order logT-\/r/logT = 
y/T logT. In d > 5, (1.7) is bounded, corroborating the classical norming. 

It is remarkable that the correlations introduced by the branching are 
strong enough that in d = 3, the limit process itself depends on the initial 
condition, not only on its density. Even though ^t, starting from 'H{'&), con- 
verges in distribution to A^, the "building up" of equilibrium is reflected 
in the different covariance structure of the renormalized occupation time 
process. 

Note that the centred Gaussian process {Xt) with covariance given by 
(1.6) can be represented as Xt = (sf^^^ + ^iY^V^, where {B[^'^^)tm is 
a fractional Brownian motion with Hurst parameter 3/4 and B^^^^^ = (see 
[3]). On the level of variances, this can be seen as follows. By (3.14) and 
(3.18), 

^ i-Nt i-Nt POO 

Var^-(Xi^)-Var^W(Xf) = ^y^ du dv J^_^^drar{0,0). 

Lemma 3.3 implies that the right-hand side is equal to Coy^^^\x!^t,Xf^) 
up to negligible terms. Hence, 

Var«('')(Xf ) p. Var^W (^(^f + X^^'. 

It remains an intriguing question to explain this representation from the 
point of view of building up a family structure in the branching random 
walk. 

For state dependent branching in d > 3, Proposition 3 in [8] shows that 

(1.8) A^ = w-lim£^W(et) 

t — >oo 

exists for any > 0. A^ is a shift-invariant equilibrium and satisfies E^'' [^0(2^)] 
i9 and E^^[{^o{x))'^] < 00. We denote a^"^ ■.= E^^[a{^o{x))], which is indepen- 
dent of X because of the shift-invariance of A^. By (1.2) and the assumption 
cr ^ 0, we have £ (0, 00). 
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In the case d = 3 with start in the Poisson field we can prove a functional 
CLT for the occupation time in zero. For d > 4, we show f.d.d. -convergence 
of the renormalized occupation time. 



Theorem 1.2 (State dependent branching). 1. In d = 3 the processes 
{X^)t>o defined in (1.4) with C{S,o) = T~i{'d) and with the norming given in 
(1.5) converge in distribution as elements o/C([0, oo),M) toward a Gaussian 
process X with covariance given by 

CoY{Xs,Xt) = p^{detQy 

(1.9) 



^1/2 cq 



^3/2 3/2_ 

2' 



^3/2 



2. In d>4 the processes (X/^)t>o defined in (1.4) with £(^0) £ {T~i{'^) , -^-d} 
and with the norming given in (1.5) converge in f.d.d. -sense toward a Brow- 
nian motion X with covariance 



CoviXs,Xt) 



(1.10) 



{27r)-^detQ)-^/^a'/ X {s At), 



2^ 



duauiO,0) + cr^'^ 



duuau{0, 0) 



X (s A t). 



d = A, 
d>5. 



Remark 1.1. For the state dependent branching case, we can at present 
not prove tightness in d > 4 due to a lack of manageable higher moment for- 
mulas. We conjecture that in the case d = 3, when starting in equilibrium, 
the limit process would again be fractional Brownian motion. In order to 
prove this using the techniques employed in the present paper, we would 
require an equivalent of the main result from [28], namely, a spatial renor- 
malization result for the equilibrium, in the state dependent case (see the 
proof of Lemma 3.2). While we believe this to be true, the techniques used 
in [28] depend on infinite divisibility, and can thus not be readily adapted. 



Corresponding functional central limit theorems for the occupation time 
of reversible interacting particle systems are well known; see, for example, 
[19, 24], or, more generally, [20] for central limit theorems for additive func- 
tionals of reversible Markov processes. In the nonreversible situation of inde- 
pendently branching system, nonfunctional versions of central limit theorems 
have been obtained in [14]. 

In order to obtain tightness in d > 4 (independent case), we have traded 
reversibility for infinite divisibility, which opens the possibility of rather ex- 
plicit calculations. 4th moment calculations are feasible, although cumber- 
some, because of the independence of families founded by different particles. 
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A program similar to ours has been carried out by Bojdecki, Gorostiza 
and Talarczyk in [6] and [7] in the case of independent branching in a some- 
what different scenario with completely different techniques: They consider 
critical binary branching particles in W^, where the individual particle moves 
according to a symmetric a-stable process, with a € (0,2], and obtain the 
following results: for a < d <2a, starting from a homogeneous Poisson pro- 
cess, the occupation time requires a nonclassical norming and converges to 
sub-fractional Brownian motion [the centered Gaussian process with covari- 
ance given by (1.6)], whereas the limit process is Brownian for d>2a (with 
a logarithmic correction to the norming in the boundary case d = 2a). Bo- 
jdecki, Gorostiza and Talarczyk have also considered the scenario with a 
heavy-tailed offspring distribution, see [4] and [5]. 

Our set-up is different in the following respect: we allow for a state de- 
pendent branching rate, we consider the lattice instead of continuous space, 
and we focus on the occupation time of a single point, whereas Bojdecki, 
Gorostiza and Talarczyk consider the integral of the particle system against 
test functions from 5(R'^). As to the techniques, Bojdecki, Gorostiza and 
Talarczyk rely on computations of Laplace functionals and Fourier analysis, 
while, in our case, the discreteness of space allows to use martingale decom- 
positions of the occupation time, and to employ techniques from the field of 
interacting particle systems (similar to [19] and [24]). Our scenario, namely, 
individual motion with a finite second moment, combined with critical bi- 
nary branching, corresponds to the case a = 2. This invites to conjecture that 
if we used an individual motion which is in the domain of attraction of an 
a-stable law (with general a G (0, 2]), we would find the same a-dependence 
of regimes as Bojdecki, Gorostiza and Talarczyk. On the other hand, our 
Theorem 1.1, part 1 suggests that in the scenario of [6], starting off from an 
extremal equilibrium for the branching system instead of a Poisson process, 
the limit process should be a fractional Brownian motion. This has, in fact, 
been proved by Milos, [23], using techniques similar to [6]. 

The rest of this paper is organized as follows: We collect some well-known 
facts about random walks and branching random walks in Section 2. Con- 
vergence and asymptotic Gaussianity of finite-dimensional distributions is 
proved in Section 3: in the case d > 4 we decompose the occupation time 
into a martingale plus an asymptotically negligible remainder term (Sec- 
tion 3.1), in the case d = 3 we "distill" a white noise out of the fluctuations 
of the particle system and represent the occupation time as an integral with 
respect to this noise (Section 3.2). In order to prove tightness, we use mo- 
ment estimates. 
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2. Preliminaries. 

2.1. Formulas related to random walks. The underlying motion process 
has generator 

Lfix)= J2 a{x,y){f{y)-f{x)). 

The continuous time transition probabilities at{x,y) solve the backward 
equation §iat{x,y) = Lat{-,y){x). We denote the transition semigroup by 
Stf{x) ■=ilyat{x,y)f{y). Let 



g{x,y) := / at{x,y)dt 
Jo 

be the Green function and 

g\{x,y) := / e~^^at{x,y) dt 



the resolvent. We denote the Green operator by Gf{x) ■=J2y9ixiy)f{y)- 
The function x ^ g{x,0) is a solution of —Lcf) = 6o and x ^ g\{x,0) a solu- 
tion of — L(j) = 5q. Define 

ut{x,y) := I as{x,y)ds, 



the Green function of the random walk killed at time t. The function (t, x) h-> 
ut{x, 0) solves {dt — L)(j) = 5o, (po{x) = 0. The Dirichlet form of the underlying 
random walk is Y.x,y"'{x,y){4'{y) - 4'{x))'^ = 2(0, {-L)(p) for (p £ £2(Z'^)- 
Note that our assumptions on a imply the local CLT 



sup ( ( 1^ ) + 1 



at{0,x) -pt{0,x) 

(2.1) 

= o(t-(^+^-2)/2), 

where is a polynomial of degree 3k and 

x^Q-^c 



pt{0,x) = {2TTt)-'^/^{detQ)-^/^ exp 



2t 



The local CLT for discrete time random walks can be found in [10] as Corol- 
lary 22.3. Prom that one can derive the corresponding result for continuous 
time. This can be done similarly to [1], page 113, where a result on the 
Galton- Watson process is transferred from discrete time to continuous time. 
Specifically, we need the following form of the local CLT: 
(2.2) at(0,0) = (2^t)-'^/2^et(Q)-^/2 + o(t-'^/2) ast^oo. 

From this, one can conclude that ||(7(-,0)||2 = /q°° ds /q°° dtas+t{0,0), so that 
g{-,0) G ^2(^'^) in case d>5, whereas ||5'i/Ar(-, 0)111 ~ ClogA^ in case d = 4. 
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2.2. Basic results on branching random walk. A convenient choice of the 
state space for a branching random walk (as well as many other "spatially 
homogeneous" particle systems), going back to Liggett and Spitzer [22], is 



H an integer- valued measure on Z"^ : ^ j{x)fi{x) < oo 



where 7 is a strictly positive function on Z'^ satisfying ^y^ida{x,y)^{y) < 
M^{x) for some constant M > 0. Note that the dependence of jt on the 
particular choice of 7 is irrelevant for our purposes, as any random {E,{x)) ^i^^d 
satisfying sup^E[^(x)] < 00 automatically has G X) = 1 irrespective of 
7. A formal construction of the independent BRW {it)t>o as an j£-valued 
Markov process can be found, for example, in Section 1 of [17]. The BRW 
with state dependent branching rate is constructed in Section 2.2 of [8]. The 
generator is given by 



(2.3) 



+ E ^^(^(r'+)+nr'-)-2F(e)) 



with ^^^y = ^ — ^x + ^y, = C + ^x a-iid = ^ — 5x- The branching random 
walk with state dependent branching rate {S,t)t>Q with initial condition ^0 £ 
X can be constructed as the unique solution to 



Ux) = eo(x) + E [ / 1(6- (y) > n)Ny^-{dsdn) 

liis^x) >n)N''^y{dsdn) 



y^tx 



[0,i]xN 

+ / 1 (is- {x) > n, ^^^} > u) N'-^'+ids dn du) 

/[0,i]xNx[0,l] V C2is-{X) J 

- [ (x) > n, > u] iV^'- {ds dn du) 

J[0,t]xNx[0,l] V C2is-[X) J 

for all X G Z"', t > 0. Here, iV^''*', x / y, are independent Poisson processes 
on [0,00) X N and N^'~^ , , x G Z*^, are independent Poisson processes 
on [0,00) X N X [0,1], all independent of ^o- N^'^ has intensity measure 
a{x,y) dt0 di, iV^'+, N^'~ have intensity measure {02/2) dt® di® du {dt, du 
are Lebesgue measures, i is counting measure). For fixed .^O) (6) is adapted 
to the filtration generated by these Poisson processes. See, for example, [8], 
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Lemma 1 and Remark 3. Define 



'[0,t]xN 

N^'^ := [ l(^s~{x)> n, ^^^l'-^^^} > u) iV^'± {ds dn du) 

J[0,t]xNx[0,l] V C2S,s-{x) J 

(with the obvious interpretations: A^^'+ counts the number of births at x, 
N^'~ counts the number of deaths at x, N^'^ counts how many times a 
particle jumps from x to y). Thus, we can rewrite 

Ux) = Ux) + Np+ - TVf + {N^ - Npy), X G t > 0. 

Immediately from the independence properties of the driving Poisson pro- 
cesses N we get the following: 

Lemma 2.1. Assume that sup^E[^o(2;)^] < oo. The compensated pro- 
cesses 

iVf'^ = iVf'^-a(x,y) fux)ds, 
Jo 

iVf'±: = iVf'±-i fa{Ux))ds 
Jo 

are pairwise orthogonal, square integrable martingales with compensators 
given by 

{N-'y)t = a{x,y) f'i,{x)ds,{N-^+)t = {N^nt = h tcT{Ux))ds. 
Jo Jo 

For /tG^i(Z^), put 
(2.4) Ft(e):=(/t,e-^A)= ftixmx)-^). 

Note that this sum is well defined if sup2,E|^(x) — 'd\ < oo. 

By compensating the driving Poisson processes, we obtain the following: 

Lemma 2.2. Let /:[0,oo) xZ'^^R satisfy supj<2-(||/i ||i + ||i9t/(||i + 
||^/i||i) < c«, and let Ft be defined by (2A). Then we have 

{dt + C)Ft{i) = {{dt + L)fui-M), 

for t £ [0, T], and the martingale Mt := Ft{Ct) - ^0(^0) - /o (^s + C)Fs{Cs) ds, 
<t<T, can be represented as 



Mt=Y f fs{x){dN^'+ - dN!'-) + E f\fs{y)-fs{x)) 



dN^'y. 
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We state some basic properties of a critical (finite variance) branching 
random walk in d > 3 which we will need in the following. Let a{x,y) = 
^{a{x, y) + a{y, x)) be the symmetrized transition kernel. (In our case a = a.) 

Proposition 2.3. Assume that a is transient. Then for each '& > 0, 
there exists an extremal invariant probability measure G V{N{'L'^)) with 
J ^(0)A^((i^) = •&. Each is translation invariant. If JZ{^q) £ {Ti.{i}), A^} , 

J t fiCs)ds -^ f fiOAm) inLi 
t Jo J 

for linear bounded, local functions f . 

Remark 2.1. In case of independent branching A^ is unique and the 
convergence also holds almost surely. 

For the independent branching case, the earliest version, in a discrete 
generations setting, appeared in [21], Satz 5.4 and Satz 5.5. A corresponding 
statement for a continuous-time model (which differs from our definition of 
a branching random walk only insofar as birth and motion are coupled) is 
contained in Theorem 6.3 and Theorem 6.4 of [15]. A proof of Proposition 
2.3 can also be obtained from the proof of Theorem 2(a), Case 1 in [17] by 
specializing to p = 0. Results in this spirit are well known, see, for example, 
the references given for Theorem in [2], which states the corresponding 
results for two "continuous relatives" of branching random walk, namely, 
branching Brownian motion and the Dawson-Watanabe superprocess. 

For Proposition 2.3 in the state dependent branching case, we refer to 
Proposition 3 in [8], where A^ is constructed as the weak limit of the system 
started in the Poisson field, and Theorem 1 in [9], which proves convergence 
of time averages. 

3. Finite-dimensional distributions. 

Proposition 3.1. Let £(^o) e {'H{'d),A^]. As N ^oo, the processes 
defined in (I. A) converge in finite dimensional distributions to a Gaus- 
sian process X (whose covariance structure depends on d, and the choice 
of the initial condition, as specified in Theorem 1.1 and in Theorem 1.2). 

The rest of this section is devoted to the proof of Proposition 3.1 in the 
various cases. 
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3.1. The case d > 4. Our strategy is as follows: similarly to the technique 
applied in [24], we are looking for a function G{(,) that satisfies LG{^) = 
(C(0) — + "small error" in order to obtain a representation of the form 

centered occupation time = martingale + "small error term" . 

We then use a general functional central limit theorem to treat the mar- 
tingale term, while we use second moment estimates to show that the error 
term becomes small. Put 

Ga(0= E 5A(x,0)(e(x)-^?), 

where g\ is the resolvent of the underlying random walk. By Lemma 2.2, we 
have 

(3.1) (Aid -£)Ga (0=^(0)-^. 

Again by Lemma 2.2, 



ds 



(3.2) := Ga(6) - GA(eo) - f CGx{is 

Jo 

(3.3) = Yl i9xiy,0)-gxix,0))Nr+Y.9xixMN^'^-Nr) 
is a martingale. Using (3.1), we obtain a representation 

[\U0) -^)ds = -Ga(6) + GxiCo) + A r Ga(6) ds + 
Jo Jo 



(3.4) 



--■.R} + Mt 



We choose A = 1/A^ and we study the terms h(i{N)~^R^^ and hd{N)^^Mpl^ 
separately in two steps. 

Martingale part: Using Lemma 2.1, we have 

(MV^)t= Yl {gi/N{y,0)-9i/N{x,0)f tci{x,y)Ux)ds 



(3.5) 

+ E 9i/Nix,0f / a{Cs{x))ds. 

Case 1: {d > 4) The martingale N~^/'^M^J^ has globally bounded jumps 
{g\{x^Q) < 5(a;,0) < ||g||oo < oo), furthermore, the jump size tends to as 
N — > oo. (3.5) yields, for any fixed t > 0, 

, 1 i-Nt 

{N-^/Hili!')t = E «(^'y)(5i/iv(y,o)-9i/^(x,o))'-y^ Ux)ds 
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1 Z"^* 



Nt 

ds 



const • t. 

This can be seen as follows. We decompose in both terms the sum over 
X € Z'^ into the sum over a ball with a large but fixed radius and the 
sum over the complement of this ball. For each point x inside the ball, 
we use that {1/N) J^^ £^s{x) ds and (1/A^) J^^ a{£,s{x)) ds converge to M re- 
spectively a^t in Li by Proposition 2.3. For x outside the ball, we estimate 
^['^(Cs(3^))] ^ C2i? and we use that g is in £2(2'^) for (i>4. Then the sum over 
the complement of the ball is small if the radius is large enough. This proves 
that the r.h.s. converges in L^, so, in particular, it converges in probability. 

Using Proposition II. 1 in [25], we conclude that {N~^/'^ M^Jj^ )t>Q con- 
verges in distribution to the law of a Brownian motion. 

Case 2: {d = 4) Here we have to slightly modify our approach because the 
Green's function is no longer in £2(Z^). Instead, we note that 

Yl 9i/N{x,^f const. > 



and that 



logiV 



a(a;,y)(5l/Ar(y,0) -5ri/7v(x,0))^ 



X 



logiV 



(ffi/Jv(-,0),(-L)9i/^(-,0)) 



2 / 1 

5i/7v(-,0),(5o - ^9i/Ar(-,0) 



logiVX""'"' ' N 
< r\9i/NiO,0) < y\g{0,0) 0. 

log A/ ' logA* N^oo 

We then argue analogously to the case above that ( . ^ M\i!^)t>n con- 

>=> J ^ ^N\ogN ^* ' - 

verges in distribution to a Brownian motion. 

Error part: Let us first consider as the initial condition. We estimate 
jgA^ [(Ga(Co))^] in order to treat the remainder term. By Lemma 3.3, we have 



^ 5A(x,0)5A(y,0)Cov^^''(eo(x),eo(y)) 

f-oo rco 

dte~^* dse"^" y at{x,0)as{x,0) 
Jo Jo 
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(3.6) 



+ ^ I dte 



ds e 



-Xs 



du ^ at{x,Q)as{y,Q)au{x,y) 

eq 



:/ dt dse-^^'^'U'dat+s{0,0) + ^ duat+s+u{0,0) 
Jo Jo I ^ Jo 



dre"^V|i9a^(0,0) + ^ j dva^{0,0) 
For d > 4, we estimate, using (2.2) 

IE^^[(GA(eo))^] <C{1+ / dre-^V(r-'^/2 ^r-'^/2+i) 



eq 



to find that 



<2C 1 + 



C 



E^-[(iV-i/2Gi/^(eo))2] < ^ + C / e-^/^r-'^/2+ 



,dr C 



■ + 



0. 



C 



The case d = 4 can be treated analogously: 

E^-[((iVlogiV)-V2Gi (^,))2]<^_ , 

' A* log A* log A* J I 

Thus, the second term of Rt in (3.4) converges to in after norming 
with hd{N), so, in particular, it converges to in probability. By the time- 
stationarity of (S^t) started from A^, we see that also the normed first term 
in (3.4) converges to in probability. Finally, the remaining integral term 
can be estimated in the following way: 

1 1 r^^ 

Gi/N{Cs)ds <tE^^|/id(Ar)-^Gi/jv(eo)| 0. 



hdiN) N Jo 



N^oo 



Putting things together, we conclude that {hd{N)~^ R^^^ )t>o — > as A^ ^ oo 
in the sense of finite-dimensional distributions. 

Now consider Poisson initial conditions. Note that 

E''(''^[(Ga(^o))']=^ E 9x{x,0f = ^ / ds / cine-"(^+")a,+„(0,0), 

a term which already appeared in (3.6). For the first term in (3.4), note that, 
by Lemma 3.3, 

= 9A(x,0)gA(y,0)Cov^('')(6(x)6(y)) 



14 



M. BIRKNER AND I. ZAHLE 



Jo Jo 

/•oo /"OO rt 

+ ds dn / dre-^('^+")a,+„+2r(0,0)E^W[a(et_,(0))]. 
JO Jo Jo 

Estimating E'^('') [a(Ct-r(0))] < C2i?, we obtain again a term which aheady 
appeared in (3.6). Since this estimate is uniform in t, we get immediately 

convergence of the last term in the definition of R^J}^ in (3.4) to zero. 

3.2. The case d = 2>. The decomposition (3.4) of the occupation time in 
a martingale term and a remainder term as for the case d > 3 can not be 
used in the case d = 2>: First, iV~^^^Gi/7v(C) does not become small in L^, 
second, as the limit process cannot be a Brownian motion, the Rebolledo- 
type arguments we used above would not help anyway. 

Our approach, again inspired by [24], is to instead "distill" a white noise 
out of the space-time fluctuations of the ergodic branching random walk 
system, and to express the normalized occupation time process as a linear 
functional of this approximate white noise. Technically, for a (momentarily 
fixed) time horizon T, we decompose the occupation time in a term Mj, and 
a remainder term, where Mj, is the final value of a martingale {Mf)t<T- 

Recall ut{x,0) = Jq as{x,0) ds and define 

Now 

(3.7) := Uf{it) - U^iCo) - [\ds + C)U^{Cs) ds 

Jo 

is a martingale, and as {dt — L)ut{-,0) = 6{-,0), we obtain, using Lemma 2.2, 
the following decomposition of the occupation time: 

riUO)-^)ds = M^ + U^{^o). 
Jo 

Being interested in iV"^/^ /o^^(6(0) - '&)ds, we find ourselves obliged to 
study N-^I^M^^ and N-^I^U^^ [ia). 

Lemma 3.2. Let C{^q) G {A^,'H(i9)} for the case of independent branch- 
ing or £(1^0) = '^{'^) fof the case of state dependent branching. The processes 

{N-^/'M^^)t>o and {N^'/'no"^^ mT>o 

converge jointly in the sense of finite- dimensional distributions to indepen- 
dent Gaussian limits. 
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Proof. We first consider {N~^/^Uq^'^ {£,o))t>o- If we start from a Pois- 
son field, that is, vC(Co) = 'H{'d), N~^/^U^'^ will converge in finite-dimensional 
distributions to the zero process: The norming with N~^/^ is too strong in 
this case, as can be seen, for example, from 

iV-3/2jE>^W[(C/o^^(eo))'] =iV-3/2^ E W(x,0)2 = 0(iV-i). 

On the other hand, if /3(Co) = the norming will be adequate. In the 
independent branching case the processes N~^/^U^'^ will have a nontrivial 
Gaussian limit. Heuristically, if we could simply replace at{Q,x) by its local 
CLT analogue, we would find 

I- NT 

= N~^/*Y.Mx)-^] as{0,x)ds 
(3.8) ^N-^/^Y.Mx)-^] 

X ^'^^(27rs)-=^/2(detQ)-V2exp(^-^^^) ds 

where ipT{x) = /(f (27rr)~^/^(det exp(- ''^f ~^^ ) dr. If, furthermore, ipr 

were a Schwartz function, we could conclude using Theorem 1 in [28]. The 
method of proof used there can be adapted to our situation; technical details 
can be found in [11], Lemma B.l. 

Now let us consider M^j. . Using Lemma 2.2, we can write [we abbreviate 
Us{x) ■.= Us{x,{))] 



^t= E f\uT~siy)-UT-s{x))dN^ 

+ J2 fuT-s{x)dN^^+- J2 fuT-s{x)dN^^-. 

Now we replace t and T by NT and multiply by N~^/^ which yields 
N-^'^M^^ = Zi{N,T) + Z2iN,T) - Z^{N,T), 

where 

rNT 

Z^{N,T)=N'''/^ J2 / {uNT-siy)-UNT-six))dN^'y, 



16 M. BIRKNER AND I. ZAHLE 

rNT 



Z2iN,T) = N'-'/^J2 / nNT-s{x)dN^'+, 

/•NT 

Zs{N,T) = N-'/^Y. / UNT-six)dN^^-. 

We proceed in two steps. In the first step we investigate Zi{N,T) and in 
the second step we consider Z2{N,T) and Z^{N,T). 

Step 1: The term Zi{N,T) converges to zero in probabihty, since the 
second moment converges to zero: 

rNT 

E[(Zi(iV,r))2]=^9iV-3/2 V a{x,y) (uNT^siv) - UNT^s{x)f ds 

pNT 

= ^Ar-3/2 a{x,y) (n,(y) - n,(x))' 

Jo 



NT 



NT 



NT 



< 2?9iV^3/2 / Us{0) ds < -dN-^/^TgiO, 0) ^ 0. 



Step 2: Now we consider the remaining terms Z2{N,T) and Z^{N,T). We 
define a random field Yat.t on L^([0,T] x M'^) via 



(Y^^T,^) ■■=N^'^ I dz r dN^^'^«'^^{s,z) 

JR^ Jo 



"N 

rT 



where lz\]\r is determined by Iz\n G 1? l\fN and z G [zJat + OAr, with f^Tv = 

Z2(iv,T) = (y^,T,w)> 

where 

Next we wish to show that Y/v^t converges toward a white noise 1^ on 
[0,r] X (with covariance measure given by times the Lebesgue 
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measure). Furthermore, for large N, the CLT suggests that vn,t should be 
similar to 

i-T-s 

vt{s,z) := / priz,0)dr, 
Jo 

where Pr{x,y) := (27rr)-3/2(det Q)-V2 exp(- iy-^)^Q-Hy-^) y Thus, we ex- 
pect Z2{N,T) ~ {Yt,vt), which shows the Gaussian nature. We proceed in 
two parts to justify these heuristics: 

Part 1: Here we show that {Yn,t^'-p) (Xt^'-p) as — > od when ip G 
L^^^2)o-°i ([0' ^ ^'^)- The index ^cr'^ indicates that this is the L^-space 

corresponding to ^cr'^ times the Lebesgue measure on [0,T] x M^. We write 
\\(p\\2 for the norm of in this space. 1^ is a space-time white noise based on 
^cr'^ times the Lebesgue measure on [0, T] x R^, that is, a. Yr :f {f, Yt) is 
a linear isometry from L'^-^^^s^^cq{[0,T] x M^) to the space of Gaussian random 

variables equipped with the L^-norm. See, for example. Chapter 1 of [27] for 
background on white noises. 

First we consider test functions consisting only of finitely many steps: Let 

n n m{k) 

(3.9) ^{s,x) = pWr fc j,,(s)1aJx) = lAk{x) l[^fc^tfc](s), 

fe=i ^'=1 ' ' ' k=i 1=1 

where yli, . . . , ^4^ C are disjoint (say, bounded closed parallelepipeds) and 
rf<tf<rt<t|<...<r^(,)<t(,).Let 

zf'^ = ivV4 ^ A(A^n(x + Ojv))iV;f^'+, fc = l,...,n. 

Then (Z/V)o<t<T = (^r\--->^t '")o<i<T is an M"-valued martingale. The 
assumptions of Proposition ILl in [25] are fulfilled since: 
(i) We observe for k ^ I that 

^zN,k^ zN,i^^ ^ ^1/2 ^ ^(^^ ^ ^ n^))A(A n (x + n^v)) 

i-Nt 

X / laiU^x))ds, 
Jo 

which converges in probability to since 

X #{x G 1? l\fN : Afe n (x + Vtj^) 7^ 0, A n (x + Hat) 7^ 0}, 
which is if is large enough since Ay^ and Ai are closed and disjoint. 
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For k = L we calculate 



pNt 



which can be seen by the following argument: 

cNt 



^ {X{Aknix + nN))f I a{U^x))ds-al''XiAk)t 



E 



(3.10) = J2 ^« / 



Aft /-Aft 



where tat— >0 as N ^ oo. Proposition 1 of [9] yields 
E'^^'''^[fiCs)giCs+uiz + ■))]- J fdA^ J gdkt 



sup 



as s, li — > oo. 



From this and shift-invariance of we can conclude that the expression on 
the r.h.s. of (3.10) converges to zero. 

(ii) We observe that Z^'^ has jumps of size A^~^/^, such that condition 
(ii) of Proposition II. 1 in [25] is obviously fulfilled. 

By Proposition II. 1 in [25], we can conclude 

(3.11) (Z^'\ . . . , Z^'") -^{Z\..., Z'^), 

A''— >oo 

where {Z^} are independent Brownian motions with variance parameter 
Icy'^^KA,). 

For (/? defined in (3.9), we obtain 

n rn(k) n rn(k) 

{yN,T, v^) = E E ^ - ^.f ) — ^ E E (4 - 

k=i 1=1 ' ' '^^°^k=i 1=1 

The limit is a sum of independent normal random variables by (3.11). There- 
fore, the limit is normal with variance J2k=i ^<^^^^(^/c)(if — and, 
hence, 

{YN,T,f) {Yt,'p)- 
Then we can extend the convergence statement to all ip G L'^^^^^ccidO, T] x 



since the functions of the form (3.9) are dense in L^-^^^ycq{[0,T] 



X 
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and since (for 99 : ^ M such that 99^ is Riemann integrable) 

lim K[{YN,T,^f] = ^crl'^ [ [ ifis, zf ds dz = \\ip\\l. 

The last assertion can be seen by the following calculation (note that 
E,[a{^j\[s{VNx))] does not depend on x): 

= Ari/2 [ ds f dzip{s,z) iViE[a(CiVs(ViVx))] 



^1/2 r ds\ [ dzifis,: 

Jo Ux+n^ 



E 



ds 



dzip{s, z) 

1 eq rT f „2 



N^a^ + TN (with TTV ^ 0) 

dw;(/?(s,i(j) 



3/2 



which is a Riemannian sum for 2'^,? /o /rs 'z' (•^j -2) c^-z ds. This completes the 
proof of the assertion. 

Part 2: Now we show that {Y^ t,vn T — vt) — > 0, in fact we will show 

' ' N—*oo 

that 

(3.12) n{YN,T, VN,T - VT?] 0. 

We have 

E[{Yn,T,Vn,T - VT?] 

= JV3/2 I dsi j dz{vN,T{s,z)-VT{s,z)) 

X ^E[a{^Nsi^x))] 



X E 



V JO 



(iz (A^^/2QAr,.(\/iVx,0) -priz,0)) 



T 

V JO 



= ^C2^N^/^ I ds[ I dr...+ I dr . . . 







2 


■ rT-s 






/ dr . . . 


+ 








Jo 
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Now note that 

^3/2 r ds\ r dr [ dzN^/'^aNri^x,0) 

r 1 r^^ 1^ 

= N^/^T J2 drar{VNx,0) 

= N~'/^T\\uNe{-,0)\\l<CV^, 
where we use for the last estimate that 

Il'"i(-)0)ll2 = y] / / dsar{x,0)as{x,0) 
xez'^-^o Jo 

= dr J ar+s(0, 0) ~ const. X \/i 

by (2.2). Similarly, 

^3/2 r ds\ r dr f dzpr{z,0)\ 

= N^/'^T V \ r dr f dzpr{z,0)\ 



<N^/^T y |f]jv| / dz( r drpr{z,0) 

= T / / dr f dspr{z,0)ps{z,0) 
Jr3 jo jo 

= 2T r dr r dspr+s{0,0) < Cy/e, 

Jo Jr 

where we used the Cauchy-Schwarz inequality. 

In order to treat the remaining term we use that [see, e.g. (2.1)] 



(3.13) 



\N'^/\Nr{VNx,0)-pr{z,0)\<Cej 



+ \z\'^/r 



uniformly in A^, r G x G Z^/\/iV, z x + Oat, where il^{N) — > as 

N ^ oo. (Note that this requires only a second moment assumption on a.) 
This yields 



x&Z3/^ 



ds 



■T-s 



dr [ dz{N^/'^aNriVNx,0) -pr{z,0)) 



.J e 



x+Qf 
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Ce 



KN'I^NfT Y: (/ drf dz- 



= C^T^i^{Nf f dz(l + |z|VT)-2 0. 



Combining, we see that limsup;v„^oQ E,[{Y]\f^T,VN,T — ^t)'^] 1^ C\/£) now let 
e ^0 to obtain (3.12). 

Thus, we have shown that Z2{N, T) converges to a Gaussian hmit. Z^{N, T) 
can be treated completely analogously, and as it involves only integrals with 
respect to {N^'~), x £ Z^, and the martingales A^^''~ and N^'~^ are all pair- 
wise orthogonal, we see that Z2{N,T) and Z3{N,T) converge jointly to (in- 
dependent) Gaussian processes. Thus, {N~'^/^M^j,) converges as ^ oo 
to a Gaussian process. 

Finally, a remark on the joint convergence of U^'^ and M^'^ when start- 
ing from the invariant distribution is in or der: Note that C/jvf (^o) de- 
pends only on the initial condition, whereas M^'^ is a function of the driving 
martingales iV^'^, x E Z^. Scrutinizing the proof the reader will find that 
even conditional on = ^) will converge to the same Gaussian pro- 

cess, as long as rj is such that >C(^t|Co = ?/) =^ as i — > oo. Note that A^-a.a. 
initial conditions have this property because A^ is an extremal equilib- 
rium in the classical branching case. The argument can be made precise 
by considering a joint characteristic functional of and U^'^ and then 

conditioning on ^o- D 

3.3. Covariance computation. First we need a second moment formula 
for branching random walks. 

Lemma 3.3. For u<v, x,y £ l/, we have the following moment formu- 
las: 

+ r a,.u+2r (x, y)E^ W [aiCuMm dr, 







K^''Mx)Uy)]=^^ + ^a^-uix,y) + ^af / ar{x,y)dr. 
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Proof. In [8], Lemma 4 one can find the moment formula, but only for 
the process in different space points at the same time. To obtain the formula 
for different times u <v, simply condition on the configuration then use 
the Markov property and the first moment formula from Lemma 4 in [8], 
finally use the second moment formula from the same lemma. Moreover, we 
use that 'E'^^^^[a{£,u~r{x))] does not depend on x. □ 

Now we compute the covariance of the limit of the renormalized occupa- 
tion time. 



Proposition 3.4. The variance of the limit of the renormalized occu- 
pation time is 



[ V2 

37r3/2 



[delQY^/^a'^it^'^ + s^'^ - \t- 3]^/% 



2^/2 



TV^oo 



d = 3,fi = A^, 

3,3/2 (det Qr'/'^T [t'^' + sf^' - lit + sf/' 

d = 3,fi = n{'d), 
{27ry^det qy^/^af X {s At), d = A, fi e {A^,n{^)}, 



roo roo 

2'd duau{0,0) +a'^'^ dntia„ (0,0) 
Jo Jo 



X (s At), 



d>5,fie{A^,ni^)}. 



Proof. The proof is split up into different cases. We assume s < t 
throughout. Let us first consider the situation C{(,o) = 7i{'d). By Lemma 
3.3, we have 



1 



(3.14) 



hd[NY Jo Jo 



Ns 



Nt 



du (it;Cov^W(Cn(0),e.(0)) 



hd{NY Jo Jo 



Nt 



du / dv a^^_u\{0,0) 



+ 



1 



Ns j-Nt puAv 

du / dv 



dr 



hd{Nf Jo Jo 

X a|„_„|+2r (0, 0)E^('') [<T{iuAv-r (0))] 

= :h+h- 

Case 1: Let d = 3. We have < Ii < dN~^''^{Ns) a^(0, 0) dr = 0{N~^/^), 
so that this term is asymptotically negligible. Fix e > for the moment. By 
(2.2), we can find If > such that a, (0,0) < (l + e)c3r-3/2 and [c7(Cr-(0))] < 
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(1 + for r>K, where cg = (27r)^3/2(jg^.g)-i/2^ rj.^^^^ bound 

h by 



2(i + e)v;;^c3 f^^-^ 



(3.15) 



iV3/2 

, (1 + erases 



du 



Ar3/2 

The first term in (3.15) is equal to 
2(l + e)V;;^C3 /-^-^ 



Ns 

du 



dv 



dr 



Nt 



{v -u + 2r)3/2 
dr 



dv 



Ns 



(?; - li + 2r)3/2 



iV3/2 

4(l+e)V^'^C3 



Ar3/2 



du 
2 



u+K 



{Ns - u) 



{v-uyl'^ (i; + u)i/2 



3/2 _ ^-1/2 



1 



f(Ar, + ^)3/2 + ^(2^ + ^)3/2 



u=0 



N—^oo o 

Analogously, the second term in (3.15) is equal to 

2(i + e)v;;^c3 f^^-K 



Ar3/2 



2(l + £)^a°qc3 

iV3/2 



d^{[(^_^)l/2_(^ + ^)l/2].=7Vt| 



.hm-uf/' + ^iNs-uf/' 



l{Nt + uf/^ + -{Ns + uf/^ 
3 3 



u=Ns-K 
u=0 



-^1(1 + efal%s{2t'/^ -{t- s)3/2 - (t + s)3/2 - (2 - 2^/^)s^/^). 

TV^oo 3 



^3/2^ 



Combining these terms and letting e — > 0, we see that 

limsup/2 < Icsa^^t^/^ + s'/^ - i(t - - ^(t + s) 

liminfTv^oo ^2 can be analogously bounded from below, concluding the proof 
in this case. 

Case 2: Let d = 4. We have < Ii < ^{NlogN)-\Ns) a^(0,0) dr = 
0(1/ log A^), so that this term is again asymptotically negligible. Arguing as 
in case 1, we can now bound I2 from above by 



2(l+e)V°% f^^-^ 



A log A 



du 



dv 



dr 



u+K Jo iv-u + 2r)2 
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(3.10) +(l±£)!f£ir-%„r.„r * 



A^logiV Jo Jns Jo {v-u + 2ry 

+ 0((logiV)-i), 

where C4 = (27r)~^(det Q)"^/^. The first term in (3.16) is equal to 
' du av ' 



NlogN Jo Ju+K \v-u v + u 



A^loffiV 



du{log{Ns — u) 

- log K - log(iVs + u) + log(2n + K)) 



= ^^^-^^;^^^i^[-{Ns - u) log{Ns -u)-u-ulogK 

- {Ns + u) log(iVs + u) 

+ u + {u + K/2) log(2n + K) - u]l=^'-^ 

— > (1 + e)^(T^'^C4S as iV — > CO. 

Now the second term in (3.16) is bounded above by 

(l + e)2a,7c4 f^'-^ , f^' dv 
' du ' 



2NlogN Jo Jns v — u 

+ {Ns - u) ]og{Ns - u) + ujl'^'-" 

= 2jVlogT' ^^^ " - ^) + ^) l°g(^(* -s) + K) 

+ KlogK - Nslog{Ns)} 

— > 0. 

Thus, letting e — > 0, we see that limsup^y^oo ^2 < cr^c^s in this case. Again, 
liminfAT^oo h can be bounded analogously, completing this part of the proof. 
Case 3: Let d>5. We have 

i-Ns ^ rNs rNt 



rNs 



2'd 
''=N. 


rNs 
'0 


2-d 


i-Ns 




'0 


= 2-ds 


Jo 



drar{0,0) +0{l/N). 



FUNCTIONAL CLT FOR BRANCHING RANDOM WALK 



25 



We decompose I2 as 

2 rNs rNs ru 

(3.17) "^^^ 

+ -/ dn / (it;/ dray^u+2r{0,0)E^^^)[c7{CuM0))]- 



The second term in (3.17) is bounded by 

(J rNt POO (J rNs-l rNt roo 

— dv drar{0,0) + — du dv drar{0,0) 

JV Jns Jv-Ns ^ Jo JNs Jv-u 

0(iV-^) + 0(iV2-d/2) 



by (2.2). 

Choose K large enough such that E'^('^)[<t(^^(0))] < (1 + e)cr.^'^ for all r > 
K. The first term in (3.17) can be bounded by 

2(l + e)cj^;;'^ r^' r^' r^-K 



N 



du dv dra^,-u+2ri0,0) + O{N-^) 

K Ju Jo 



= {l + e)aT du dv drar{0,0) + 0{N-^) 

Jk/n Jo Jv 

/•oo /"OO roo 

— > {I + e)a'^s / dv dr ar{0,0) = {1 + e)a^s / duuau{0,0). 
Jo Jv Jo 

Letting e — > 0, we see that limsupjv^^oo of the first term in (3.17) is bounded 
by a^s /q°° duuau{0,0). Since the liminf can be bounded analogously, we 
are done with the case d > 5. 

Now let us consider the situation -C(^o) = Aij. Then we have, by Lemma 
3.3, 

,Q pNs pNt 

jgA, [x^x^] = —— j duj dv a|,_„| (0, 0) 



^cq /.Ns rNt 



(3.18) + ^, / du dv drariO,0) 

2hd[Ny Jo Jo J\v-u\ 

= :h + l2- 

The computations for d = 4 and d>5 are entirely analogous to those above, 
and will be omitted. Let us briefly comment on the case d = 3 in this situa- 
tion. Ii is again negligible, and choosing K large enough, we can now bound 
1 2 from above by 

2(l + e)a;;^C3 f^'-'' dv 



iV3/2 7o 7^+^ _ ^)l/2 



26 



M. BIRKNER AND I. ZAHLE 



+ 




We conclude the proof as above. □ 

4. Proof of Theorem 1.1. Here we complete the proof of Theorem 1.1. In 
view of Proposition 3.1, it suffices to check that the sequence , G N, is 
tight (e.g., in the space of all continuous processes, equipped with the norm 
of locally uniform convergence). In order to do so, we use the well-known 
criterion on moments of increments, stating that a sequence of processes 
is tight (and, furthermore, any limit point has continuous paths) if 
there exist a,f3 >0 such that, for each to > 0, 



holds uniformly in N and < s < t < to (see, e.g., [18], Corollary 14.9). 

For d > 4 and independent branching, with a = 4 and P = 1, this is the 
content of the following lemma. 

Lemma 4.1. Let d> 4 and ^ € {Ti.{'d),A^} . Then in the independent 
branching case for each to > there is a C = C{to,'d,d) such that, for all 
< s < t < to, we have 



Obviously, this lemma requires 4th moment computations. We do this by 
a "bookkeeping of trees" similarly as in [16] for superprocesses. We refer to 
[11], Lemma 4.1 for details. See also [7], proof of formula (3.36) for a related 
approach using Laplace transforms. 

In the case d = 3 it turns out that second moments [a = 2, (3 = 1/2 in (4.1)] 
suffice. The corresponding estimate is provided in the following lemma. 

Lemma 4.2. Let o? = 3 and fj, G {TC{'d),A^}. For each to > 0, there exists 
a constant C = C{tQ,'&) such that 



(4.1) 



E[|Xf -Xf |°]<C7(t-s)i+'^ 




d = A, 
d>5. 



E'^[(Xf - X^f] < C{t - s)3/2 VO < s < t < to 



holds uniformly in N . 
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Proof. We begin with the independent case. Note that for any initial 
distribution fi we have 



Ar3/2 



Nt pNt 

du / d7;Cov^(Cn(0),e.(0)), 

Ns JNs 



thus, we see from Lemma 3.3 that E«('')[(Xf - X^f] < E^^[{X^ - X^)^] 
and it is, hence, sufficient to consider the stationary initial distribution A^. 
By stationarity, we can assume without loss of generality that s = 0, t < 
to. Put ip{r) := Cov^^(er(0),^o(0)). We have < ip{r) < C(l A r-^/^) by 
Lemma 3.3 and (2.2). This allows to estimate (and the same is true in the 
state dependent case) 



E^^[{X^^ - X^f] 

Nt 



du 



Nt 



dvip{v — u) 



iV3/2 

:2iV^/2 r r dv^{N{v-u)) 



Jo 

< 2CN^/'^l^t'^l{Nt < 1) + tl{Nt > 1 
. 2Ct^/^\{Nt)^/^l{Nt < 1) + {Nty'^/'^l{Nt > 1) 



1 

N 



ds—j= 

i/N VNs 



+ {t/N)-'/H{Nt > 1) (2{t/N)'/' - I) I 

For the state dependent case, we observe (again by Lemma 3.3) 
1 



1 



hd{Nf Jns 



Nt 



Nt \ 2- 

{U^)-^)du 

Ns 
Nt 



du dz;(E«W[e„(0)e.(0)]-i92 



hd{N)^ Jns 



Nt 



du 



Ns 
Nt 



u 



dv (^av~u{x,y) 

+ /" a,-u+2r{x, y)E^^^^ k(e«-r (0))] dr, 
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2 fNt pNt / nu \ 

~lhi{NYjN J dv\da^-u{x,y) + C2^ ay-u+2rix,y) drj , 
which can be estimated as in the independent case. □ 
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